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Abstract 

Plasma velocity and magnetic field measurements from the Voyager 2 mission are used to study 
solar wind turbulence in the slow solar wind at two different heliocentric distances, 5 and 29 
astronomical units, sufficiently far apart to provide information on the radial evolution of this 
turbulence. The magnetic helicity and the cross-helicity, which express the correlation between 
the plasma velocity and the magnetic field, are used to characterize the turbulence. Wave number 
spectra are computed by means of the Taylor hypothesis applied to time resolved single point 
Voyager 2 measurements. The overall picture we get is complex and difficult to interpret. A 
substantial decrease of the cross-helicity at smaller scales (over 1-3 hours of observation) with 
increasing heliocentric distance is observed. At 5 AU the only peak in the probability density of 
the normalized residual energy is negative, near -0.5. At 29 AU the probability density becomes 
doubly peaked, with a negative peak at -0.5 and a smaller peak at a positive values of about 0.7. 
A decrease of the cross-helicity for increasing heliocentric distance is observed, together with a 
reduction of the unbalance toward the magnetic energy of the energy of the fluctuations. For the 
smaller scales, at 29 AU the normalized polarization is small and positive on average (about 0.1), 
but it is zero at 5 AU. For the larger scales, the polarization is low and positive at 5 AU (average 
around 0.1) while it is negative (around - 0.15) at 29 AU. 
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1. Introduction 

The solar wind offers a unique scenario where in-situ spacecraft observations can unveil 
many aspects of very high Reynolds number magnetohydrodynamic turbulence in a magnetized 
plasma. One important issue of solar wind studies is to ascertain the dynamics of its fluctuating 
magnetic and (plasma) velocity fields. These studies involve the measurement and analysis of 
features which range from large-scale structures originating in the solar corona down to the 
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kinetic scales of the ion Larmor radius. Since the earliest observations of the solar wind mm 
it has been clear that the solar wind has strong variations in velocity, density, magnetic field and 
temperature which cannot be described only in terms of the superposition of outward propagating 
Alfven waves. The general picture (see, e.g., Q) is that solar wind fluctuations are generated at 
the solar surface, where new structures continuously appear and fade away, and then propagate 
outwards through large wavelength Alfven waves. Then energy is gradually transferred toward 
smaller scales by nonlinear processes until it goes into plasma particles in the form of heat. It is 
believed that stream-shear instabilities play a major role in the triggering of this transfer process. 

The observation of the variation of the solar wind parameters with the distance from the sun 
is essential for understanding the evolution of the interplanetary fluctuations and has been made 
possible by the launch of spacecraft. In particular, the data from the Voyager mission offer a 
good opportunity to obtain spectral information and investigate the structure of solar wind and 
its driving mechanisms (underlaying dynamics) close to the solar equatorial plane over a large 
range of heliocentric distances and solar wind conditions and also beyond the termination shock 
and towards the edge of the heliosphere. 

In this paper we characterize the dynamical structure of the solar wind magnetohydrody¬ 
namic turbulence at distances of 5 and 29 astronomical units (AU) from the Sun in terms of both 
the distribution of cross-helicity and magnetic helicity using data from the Voyager 2 mission. 
The spectral analysis of these quantities provides information on the plasma evolution over a 
large range of scales which, in the end, describe the overall macroscopic state of the solar wind. 
The magnetic helicity and the cross-helicity both are important characteristics of the solar wind 
turbulence. Like the kinetic helicity and the total energy, both quantities are invariants of the 
flow when diffusivity and resistivity are not taken into account. A non-zero magnetic helicity 
indicates the lack of reflectional symmetry in the flow and is related to the dynamo effect. The 
cross-helicity is proportional to the correlation between velocity and magnetic held fluctuations 
and measures the relative importance of Alven waves in the global fluctuation. As many studies 
in literature have shown (see, e.g. naisi), the normalized cross-helicity measures the coupling of 
the magnetic and plasma fields and thus is correlated to the self-production of turbulence. Many 
studies in literature (see, e.g. mil) indicate that the solar wind turbulence evolves by gradu¬ 
ally reducing these correlations and by showing a predominance of the energy associated to the 
magnetic fluctuations with respect to that of the velocity fluctuations. 

The definition and the significance of the magnetic helicity and cross-helicity are reviewed 
in section]^ In section]^ we present the results obtained from the Voyager 2 time series using 
reconstruction methods for spectral spectral analysis of lacunous data previously tested on I?). 

2. Helicity and cross-helicity 

Since the seminal paper by Moffatt [|8], the concept of helicity has been used as a tool to 
describe the structure of a fluid or plasma flow. From a mathematical point of view, the helicity 
associated with a vector held u in a domain V (bounded or unbounded) is defined as the integral 
of the scalar (dot) product of the vector and its curl, that is u ■ (V A u) dV. In a plasma flow, 
we can define the kinetic helicity associated with the velocity held u and its curl, the vorticity 
held a>, and the magnetic helicity //„, associated with the magnetic potential A and its curl, the 
magnetic held B: 
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The quantities ht - n ■ a> and h„, = A ■ B represent the kinetic and magnetic helicity density 
within the flow domain. Both the helicities, and //„, as well their densities, and h^, are 
pseudoscalar quantities, that is, they change sign switching from a right-handed to a left-handed 
frame of reference. Although the origin of the concept of helicity can be traced to the works 
by Kelvin 191 on the dynamics of vorticity in an inviscid flow, and by Woltjer IfTOl on ideal 
resistiveless magnetohydrodynamics, it was first explicitly introduced in more recent times by 
Moffatt in 1969 lH. Its name is derived from the fact that, for a given enstrophy |u>p, helicity is 
maximum when the vorticity oi is parallel to the velocity u, a situation in which the streamlines 
are locally right-handed (positive h^) or left-handed (negative h^) helices about the a> axis with 
a pitch proportional to hklaP'. Helicity is related to the flow kinematics since it gives a measure 
of the linkage of vortex lines of the flow. By assuming, for example, that the vorticity is zero 
in the whole domain except inside two closed and untwisted but reciprocally intertangled vortex 
tubes, it can be shown (see, e.g., ifTTlfT^fTSl ') that the helicity is equal to ±2nrir2, where Tj and 
r 2 are the circulations of the two vortex tubes (vorticity flux) and n is the number of windings 
of the two tubes (see figure [T]|. More generally, IfTTlfTTlfTSlfT^ . if we assume that the vorticity 
(or the magnetic field) is concentrated in a finite number of closed tubes, the total helicity can be 
written as the sum of the contribution of the helicities of the individual tubes, coming from their 
twisting, and of a contribution of their mutual helicity due to their winding number, a measure of 
the winding of the tubes about each other. Although, in general, a vorticity or magnetic potential 
distribution cannot be simply decomposed into a finite number of simple non-overlapping flux 
tubes, this simple sketch helps to understand the helicity kinematic significance. 

In the same way as the kinetic helicity is a material invariant of the flow in an inviscid (neu¬ 
tral) fluid, which is a direct consequence of Kelvin’s theorem, the magnetic helicity is an 
invariant of ideal inviscid and resistiveless magnetohydrodynamic flow as first proven by Wolt¬ 
jer GO). The proof of this invariance for ideal magnetohydrodynamic equations depends on the 
boundary conditions for u and B and requires vanishing surface flows on the boundaries. Under 
such conditions, magnetic helicity is one of the three quadratic invariants (often called “rugged” 
invariants), together with the total energy £ = H- Z?^/(jUop)j dV and the cross-helicity black 
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From a geometric point of view, the cross-helicity gives a measure of the degree of linkage of 
the vortex and magnetic flux tubes within the flow. In fact, if one again considers the simplest 
archetypical sketch, where vorticity and magnetic fields are concentrated in two thin untwisted 
closed tubes with are intertangled, the cross-helicity is equal to InF^, where n is the number of 
linking or winding of the two tubes, and T and are the vortex and magnetic flow across the two 
tubes. Therefore, a normalized cross-helicity gives the measure of the knottedness of the vortex 
tube(s) with the magnetic tube(s). 

While the physical meaning of energy is clear enough to omit further description, a brief 
discussion on the dynamical importance of H„, and He beyond their kinematic description can be 
helpful. Both quantities and He, as well their densities, h„, and he = u B, are pseudoscalar 
quantities, that is, they change sign under a reflection of the coordinate system. This implies 
that they are zero in a mirror-symmetric system and non-zero values can appear only when this 
symmetry is broken in the flow. Because the value of these inviscid (ideal flow) invariants cannot 
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Figure 1: Sketch of a configuration of linked vortex/magnetic tubes with different number of windings n, that produce 
a non-zero helicity. Kinetic helicity: red and blue tubes are vorticity tubes. Magnetic helicity: red and blue tubes are 
magnetic tubes. Cross-Helicity: the red tube is a vortex tube and the blue tube is a magnetic tube (or the opposite). 


be modified by nonlinear terms, they constrain the overall dynamics and, therefore, their value 
and spectral distribution can give valuable information on the dynamics of turbulence in the 
solar wind. Of particular interests are the magnetic and cross-helicities coming from plasma and 
magnetic field fluctuations. 

The magnetic helicity plays a central role in the dynamo effect, the so-called a-dynamo. In 
fact, the presence of a magnetic field without reflectional symmetry, that is, with a non-zero mag¬ 
netic helicity, generates fluctuations of the magnetic field, which, to first order approximation, 
are proportional to the helicity ( ifT^fTTl ISlll. Compared to the magnetic energy, which appears 
to be transfetTed to the small scale fluctuations, magnetic helicity presents an inverse cascade, 
that is, it cascades toward the larger scales of the flow (see mi). As a consequence, magnetic 
helicity is depleted much more slowly than energy in a freely evolving flow, while in a forced 
flow small-scale helical forcing can produce large-scale magnetic fields. 

The turbulent cross-helicity is directly related to the coupling coefficients for the mean vor¬ 
ticity in the electromotive force and for the mean magnetic-field strain in the Reynolds stress 
tensor 0. This relationship suggests that the cross-helicity is important where inhomogeneities 
in the flow and in the magnetic field appear. Since such large-scale structures are ubiquitous in 
astrophysical phenomena, cross-helicity is expected to play a role in such flows. In the presence 
of cross-helicity, the mean vortical structures contribute to the electromotive force; the genera¬ 
tion of magnetic field due to this effect is called the cross-helicity dynamo (or j0-effect). In fact, 
in presence of a large scale vortical motion, the turbulent cross-helicity contributes to the elec¬ 
tromotive force aligned with the large-scale vorticity. Provided that the velocity and magnetic 
field fluctuations are correlated, a contribution to the electromotive force parallel (when he > 0) 
or antiparallel (when he < 0) to the mean vorticity arises m. Moreover, more recently direct 
numerical simulations of magnetohydrodynamic turbulence have shown that, in association with 
high values of the cross-helicity, a blocking effect on the spectral transfer of energy is observed 
and results in energy accumulation in the system. This is concomitant with an increase of the 
exponent of the spectrum with the cross-helicity level. The spectral exponent increases toward 
the limiting value of 2 ll20l . 

It should be noted that decaying incompressible magnetohydrodynamic turbulence, in the ab¬ 
sence of large-scale mean shear, tends to amplify the cross-helicity due to a phenomenon called 
“dynamic alignment” ||2T1 which tends to increase the correlation between the velocity and mag¬ 
netic field by locally aligning the directions of their polarizations while the energy decays, an 
effect which is stronger at smaller scales. However, observations indicate that this effect does not 
play a major role in the solar wind, where a gradual reduction of the content of Alfvenic (cor- 
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1979 

DOY 1-180 

1989 

DOY 5-100 

r 

[AU] 

4.48 H- 5.28 

28.03 H- 28.94 

{u) 

[km/s] 

454 + 43 

464 + 39 

{B) 

[nT] 

0.981 

0.207 

Ek 

[km2/s2] 

1.224- 10^ 

1.005 ■ 10° 

Eh 

[km2/s2] 

O 

00 

2.341 ■ 10° 

K 

[km2/s2] 

1.58 ■ 10' 

-3.04 ■ 10° 

Ei:„ 

[km2/s2] 

4.16-10'° 

-9.79 ■ 10'° 


Table 1: Solai* wind data: r is the heliocentric distance, {u) and {B) are the mean values of the modulus of the velocity and 
the magnetic field, Ej^ and Em are the turbulent kinetic and magnetic energy per unit mass, {u'.u'pjl and {b'.b'pjl, and 
and are the mean values of the magnetic helicity and cross-helicity of the fluctuations. All statistics are computed 
by averaging the time series over the periods indicated, except the magnetic helicity which is obtained by integrating its 
spectrum. 
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o-R 



1979 

1989 

1979 

1989 

1 hour 

9.63 ■ 10 ° 

4.43 ■ 10-° 

-0.408 

-0.133 

3 hours 

O 

o 

1 

5.59 ■ 10-° 

-0.459 

-0.194 

12 hours 

1.06- 10-' 

6.09 ■ 10 ° 

-0.460 

-0.322 

24 hours 

9.67 ■ 10-° 

5.88 ■ 10-° 

-0.405 

-0.404 


Table 2: Average values of the normalized cross-helicity (Tc and of the normalized residual energy (Tr of the fluctuations 
as a function of the averaging time. 


related) fluctuations with the heliocentric distance is observed, in particular near the equatorial 
plane Il22l . 

3. Voyager 2 data analysis 

3.1. Voyager 2 data 

The present analysis uses Voyager 2 data from the solar wind plasma and magnetic field ex¬ 
periments to analyse the solar wind turbulence in the outer heliosphere Il23l . These experiments 
provide data on the plasma velocity, the density, and the magnetic field. The three components 
of the magnetic field B and of the plasma velocity u are given in the heliocentric spherical RTN 
coordinate system {eR,eT,eff} where e* points in the radial direction outward from the Sun, Cj 
is in the plane parallel to the solar equator and positive in the direction of solar rotation, and 
Cn = Cr A Cr points toward the heliographic north completing the right-handed coordinate sys¬ 
tem. The (cr, Cr) plane is inclined with respect to the equatorial plane by an angle equal to the 
current latitude of the spacecraft. For our study we have selected two periods: the first six months 
of 1979 and the first three months of 1989, as shown in table [T] The dataset consists of the ve¬ 
locity and magnetic field components and the density from 1-1-1979 00:00 GMT to 29-5-1979 
19:00 GMT for a total of about 180 days and from 5-1-1989 to 11-4-1989 for a total of 96 days. 
During 1979, the plasma velocity was sampled each 96 s and the magnetic field data are 48 s 
averages, while during 1989 the plasma data were sampled each 192 s and the magnetic field 
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Figure 2: Probability density function of the normalized cross-helicity g{(Tc) and of the normalized residual energy q{(Tr). 
(a,c) 1979 data DOY 1-180, (b,d) 1989 data, DOY 5-100. The different cui-ves refer to 1, 3, 12 and 24 hour averages. 
Panels (e) and (f) show the joint probability density function of cTc and (Tr computed from hourly averages. The circle 
indicates the upper bound cr^ + (t^^< 1. 
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data are 48 s averages. These two periods precede the Voyager 2 encounters with Jupiter and 
Neptune in July 1979 and August 1989. The Voyager 2 heliocentric distance varies from 4.48 to 
5.28 astronomical units (1979) and from 28.03 to 28.94 AU (1989) and the solar wind shows a 
good degree of homogeneity. 

During both periods the solar wind is in the so-called ’’slow wind" state, with a mean ve¬ 
locity of less than 500 km/s and a fluctuations root mean square of about 25 km/s. The slow 
wind condition is characterized by relatively high fluctuations, with a high fraction of the total 
fluctuation energy per unit mass in the magnetic energy, which accounts between 54% and 70% 
of the energy of the fluctuations. During these two periods, the mean velocity is almost aligned 
with the radial direction, while the magnetic field tends to have a strong component normal to 
the radial direction. 

We characterize the structure of the solar wind in terms of its cross-helicity and magnetic 
helicity. The determination of the kinetic and magnetic helicity, as well as the cross-helicity, 
requires knowledge of the velocity, vorticity, magnetic field and its potential in the whole flow 
domain. In the solar wind, however, the magnetic and plasma velocity fields are measurable 
only at the spacecraft position, like the Voyager mission. This restriction implies that only one¬ 
dimensional spectra can be determined by using the Taylor hypothesis of frozen flow. However, 
a major problem in the spectral analysis of the Voyager mission data is the sparsity of the data, 
which become more and more lacunous as the spacecraft moves away from the Earth. During 
the 1979 period, 28% of the data are missing for the plasma velocity and 24% for the magnetic 
field, while during the 1989 period the missing data increase to 61% for the plasma velocity and 
60% for the magnetic field. In order to determine spectral quantities it is necessary to use a data 
reconstruction method. In the present work, we have adopted data interpolation to fill the gaps 
in the time series. This method gives correct estimates of the spectral exponents with an error up 
to 5%. The error estimate was done in Gl on Voyager 2 1979 data. In this work, we use again 
Voyager 2 1979 data plus Voyager 2 data from day 5 to day 100 of 1989. We have verified that 
in this last case the same estimate is valid. 

From the observed time series we define in each time interval the mean and the fluctuating 
field. Because of the super-Alfvenic flow, we can use the Taylor hypothesis of a frozen-in flow. 
This approximation is generally accepted in Solar Wind observations (see for instance Matthaeus 
and Goldstein 1982 ll24l and Roberts and Goldstein 1987 ll25]l l. In particular, a recent paper by 
Howes, Klein and TenBarge 2014 ll2^ shows the validity of the Taylor hypothesis for super- 
Alfvenic flows. This validation is done under the premise that the frequency of the turbulent 
fluctuations is well characterized by the frequency of the linear waves supported by the Solar 
wind plasma. As a consequence, it the Taylor hypothesis holds for linear kinetic wave modes 
in the weakly collisional Solar wind plasma, with the exceptionn of high wavenumber whistler 
waves. This assumption allows us to transform time into space by the transformation x - Tit 
and, consequently, to determine spatial spectra by linking the frequency to the wavenumber, 
K = 2nf 111, where m is the mean velocity. In the following, we use the Taylor transformation by 
assuming the mean velocity field as constant in the two spatio-temporal ranges, the first six month 
of 1979 and the first 3 months of 1989. As can be seen in Table 1, the velocity field was 454 + 43 
km/sec in 1979 and 464 + 39 in 1989. The magnetic field is expressed in Alfven units, that is, in 
terms of b = B/(pop)*'^^, where p is the plasma density and pQ the magnetic permeability of free 
space. In this way, b has the dimension of a velocity and, consequently, the cross-helicity has the 
same dimensions of an energy per unit mass while the magnetic helicity has the dimension of 
an energy per unit mass divided by a length. The magnetic and cross-helicity spectra have been 
computed from single point measurement by the Voyager 2 spacecraft using data interpolation to 
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reconstruct the missing values (see El for a discussion of the errors introduced by reconstruction 
methods in the evaluation of the spectra). We applied the formulae hrst proposed by Matthaeus et 
al. 11241m to deduce the magnetic helicity and cross-helicity spectra from the Fourier transforms 
of the velocity-magnetic field correlations. For the magnetic helicity of the fluctuations 

h^{K)^^lm(STN(K)), (4) 

K 

where S jn is the Fourier transform of the correlation between the tangential and normal com¬ 
ponents of the magnetic held (S tn(k) - ^ f Crv(Oe'*'^dr, with CrNir) - {b'j{x)b'i^(x -t- r)^), 
while the cross-helicity of the fluctuations is computed as the Fourier transform of the correla¬ 
tion + u'jb'j + u'f^b'j^ between the velocity and magnetic fields. Only the assumption of ho¬ 
mogeneity is necessary but no assumption of isotropy is made in obtaining the one-dimensional 
energy and helicity spectra. 

3.2. Results 

In this study the solar wind turbulence is characterized in terms of two parameters, the mag¬ 
netic helicity and the cross-helicity /z' of the fluctuations. The last, which is defined as 
h' = (u' ■ b'), is first considered in the following. The cross-helicity can be rewritten in terms of 
the normalized cross-helicity ac 


a-c 


2h' 


<u' ■ u') + <b' ■ b'> 


(5) 


and of the residual energy ctr 


<u' ■ u') - <b' ■ b') 
<u' ■ u') -I- <b' ■ b')’ 


( 6 ) 


The normalized cross-helicity of the fluctuations, cTc, is a function of the correlation between 
the plasma velocity and magnetic held fluctuations. It becomes equal to the correlation function 
C(u',b') = (u' ■ b')/((u' ■ u')^^^(b' ■ b')*'^^) when the turbulent kinetic energy per unit mass 
(u' ■ u')/2 is equal to the magnetic energy of the fluctuation per unit mass (b' ■ b')/2, a situation 
which occurs only when the flow is isotropic. From its definition, cr^ is bounded to values 
between -1 and H-1, which are attained only if the fluctuations of velocity and magnetic field 
are perfectly correlated, a situation which corresponds to outward propagating Alfven waves 
{(Jc - 1) and inward Alfven waves (cT;, = -1). The imbalance between the kinetic and magnetic 
energy can be expressed by the Alfven ratio r ,4 = <u' ■ u')/(b' ■ b'), that is the ratio between the 
kinetic and magnetic energy, or by the so-called normalized residual energy Il271l25l . 

The absence of magnetic fluctuations corresponds to a normalized residual energy equal to 
1, the absence of kinetic fluctuations corresponds to a normalized residual energy equal to - 
1, while zero implies equipartition of the energy between kinetic and magnetic fluctuations. 
The correlation is linked to the cross-helicity through the residual energy ctr by C(u',b') = 
crc/(l - Therefore, the cross-helicity does not generally equal the cross-correlation and 

they coincide only if there is equipartition of energy between the two fields. 

Past analysis of solar wind turbulence in the inner heliosphere ^ O indicates that the solar 
wind evolves toward a decreasing correlation between velocity and magnetic held fluctuations 
and toward an increasing imbalance in favour of the magnetic contribution to the energy of the 
fluctuations. 
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Figure 3: Variation of the cross-helicity with heliocentric distance. Filled symbols show the results of the present analysis 
of Voyager 2 data (in tabular form in table|^. The lines represent the mean field models by Breech et al. (28l (continuous 
line) and by Matthaeus et al. ESI (dashed line), empty symbols are Voyager and Helios data from (23 


Present analysis shows that the greatest imbalance occurs in the less spatially evolved (closest 
to the Sun) 1979 solar wind, while the more evolved (furthest from the Sun) 1989 solar wind 
shows a tendency toward energy equipartition, at least at the small scale level. In fact, while 
the 1989 24-hour average data show a residual energy cr* of about -0.4, the hourly average 
data indicates a drop to -0.13 (see table |^. This behaviour was not present in the 1979 data, 
where all scales had a similar ctr, ranging from -0.4 to -0.46. The imbalance in favour of the 
magnetic energy tends to be less pronounced at small scales. This relative increase in magnetic 
energy could be due to waves generated by pickup ions, which initially have an unstable ring 
distribution (see ll29l l. 

Panels (a) and (b) of figure|^show the probability density function of cr^ computed using dif¬ 
ferent averaging times in (j^, from 1 hours to 24 hours. The probability densities have been de¬ 
duced from the numerical derivative of the empirical cumulative distribution functions Eomi]. 
Distributions become more symmetric with increasing heliocentric distance, with no qualitative 
difference between the different averaging times apart from a small reduction of the variance of 
cTc with increasing time scale while the mean values remain almost constant (see table |^. The 
most significant result is the strong reduction of the normalized cross-correlation variance in the 
1989 data, which indicates a marked reduction in the presence of Alfvenic fluctuations. The cr^ 
values tend to a symmetric distribution centered around zero. An increasing variance at small 
scales is observed in 1979, see panel (b) in figure]^ an increasing variance is a common char¬ 
acteristic of the inner heliosphere 1^ where broader probability density functions are usually 
observed for averages below one day. Panels (c) and (d) of figure]^ show the probability density 
function of the residual energy ctr. Positive values are almost completely absent during 1979 at 
5 AU, with a peak around -0.5 which sharpens with increasing averaging scale. On the contrary, 
the 1989 distributions at 29 AU are more scattered with a non-negligible data with positive <jr, 
in particular for the shorter averaging times, one and three hours, where a secondary peak around 
0.5 is clearly visible. Both positive and negative ctr are not correlated with a particular value 
of cTc, which remains symmetrically distributed around zero for all values of ctr (see the joint 
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probability distribution in panel (f)). 

The overall evolution of cr^ with heliocentric distance is summarized in table|^and in figure]^ 
Moving from 5 to 29 AU close to the ecliptic plane, our data show a moderate decline of ctc which 
follows well the trend predicted by previous investigations ESlI^ based on mean field Reynolds 
Averaged MHD simulations where turbulent trasport was taken into account. It should be noted 
that time variations due to the solar wind activity should not have a major influence on present 
results because both periods, 1979 DOY 1-180 and 1989 DOY 5 -100, show a slow wind and the 
time spans are longer than the solar rotation, about 25 days at the equatorial plane. In figure]^ we 
show also the normalized cross-helicity computed by Robert et al. in 1987 ||25]| (data from Helios 
1 at 0.3 AU and Voyager 2 1985 from 2 to 20 AU), see open symbols. These last computations 
agree to a lesser extent to the above mentioned Reynolds averaged simulations. However, if one 
interpolates data with a same average time interval, one sees that they follow better the older 
2004 version of the model Il22l than the updated 2005 version of the same model Il28l . In the 
last case, the transport of cross helicity and radial evolution of Alfvenicity was computed with a 
slightly lower value of the Karman-Taylor constant and slightly different boundary condition at 
0.3 AU (similarity lenght of 0.01 AU against 0.025, fluctuation energy of 2200 instead of 2000 
km^/s^ and a slightly higher normalized cross-helicity). It should be noted that in the data 
were not all properly sorted, for example by wind speed. 

As the plasma moves away from the Sun, the velocity and magnetic fields become less cor¬ 
related. This decrease in correlation is in contrast with results obtained by means of numeri¬ 
cal simulations of homogeneous and isotropic incompressible magnetohydrodynamic turbulence 
(e.g. ll^ ) and by phenomenological interpretation carried out on the same system configura¬ 
tion ( Il2n[^ ). In short, the relevance of the dynamic alignment in compressible MHD turbulent 
flows, and in particular in the the Solar wind, is still debated (see, e.g., ||3] and Il34l and the 
references therein). Compressibility, velocity shears and density gradients (and related turbulent 
diffusion) present in the Solar Wind may overwhelm the effect of the dynamic alignment, driving 
the turbulence toward lower levels of cross-helicity. 

A closer look at the scale dependence of the cross-helicity can be obtained by looking at its 
one-dimensional power spectrum. The cross-helicity spectrum can take both positive and nega¬ 
tive values. The sign of the cross-helicity spectrum often indicates the propagation of Alfvenic 
waves outwards from the sun (positive helicity) dm, as also contemplated by Belcher & Davis 
im, but periods with mixed helicity sign can appear. In our data, see figure]^ there is no preferred 
sign in any part of the cross-helicity spectra, whose modulus shows a fair power law scaling with 
an exponent equal to -1.59, close to the spectral exponent of the total energy spectrum, -1.62 
during 1979. In the later period, 1989, the spectra become steeper, and the spectral exponents 
are -1.74 for the cross-helicity and -1.80 for the total energy.This is clearly visible in figure]^ 
which shows the cross-helicity spectrum compared with the total energy spectrum, the sum of 
the kinetic and the magnetic energy of the fluctuations. The number of frequencies with positive 
and negative values is comparable at all scales, and they balance; the integral cross-helicity can 
be estimated as 15.8 km^/s^ during 1979 and 304 km^/s^ during 1989, much lower than the total 
energy of 2640 and 2346 km^/s^ in 1979 and 1989, respectively. Globally the normalized cross¬ 
correlation is less than 1% at 5 AU and of the order of 10 percent at 29 AU and therefore the two 
fluctuationg fields are very mildly correlated. 

Similar to the spectrum of the cross-helicity, the spectrum of the magnetic helicity, which has 
also a power law range with exponent -2.66 in 1979 and -2.92 in 1989, see panels (a) and (b) in 
figure]^ shows almost an equipartition positive and negative values throughout the inertial range 
of the spectrum. The exponents of the magnetic energy spectrum in the same range of frequencies 
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Figure 4: Wave numbers spectra of the cross-helicity density and of the total energy. These spectra are obtained by 
means of the Taylor hypothesis, see section]^ The mean speed is that computed over the whole time interval, see Table 
1. For claiity, in panel (a) not all the points of the spectra are shown for frequencies above 10“^ Hz. The same in panel 
(b) above 10“^ Hz. The small peaks at a frequency of 2.6 • 10“^ Hz are likely instrument-related, since such frequencies 
correspond to about four time the sampling period of the plasma velocity. 


are -1.63 and -1.80, respectively (notice that in the ordinate the quantity Kh,n is represented). This 
result implies that the relative magnetic helicity of each scale, given by the ratio a'^,„(/(')/£’„j(/c), 
slowly decreases as (5 AU) and (29 AU) in the observed wavenumber range. The 

X +oo 

hm{K*)6.K* carries the information on the 
polarization at the different scales. For instance, the sign of first the points of the dark curves 
in panels in panels (c) and (d) gives the global polarization because of course the largest scales 
contribution prevails energetically. The sign of H,„l {E^Ik*)Ak* is negative at the largest 
scales during 1989 but is positive during 1979. If we consider all the wave numbers, this value is 
relatively small: 0.135 at 5 AU and -0.132 at 29 AU. But if we consider the small scale ranges, 
we see that in 1979 the polarization is negligible, while in 1989 it is definitely positive for the 
last two decades. It is known that magnetic helicity can change sign and have opposite signs at 
the largest and smallest scales in magnetohydrodynamic system, an effect of the inverse cascade 
and its conservation 1351 . As a general comment, the weak values of the normalized magnetic 
helicity can be considered a signature of the negligible injection of helicity in the solar wind at 
the equatorial plane because the presence of the inverse cascade of helicity makes the helicity 
dissipation much lower than the dissipation of energy, thus the magnetic helicity remain almost 
constant throughout the system. 

4. Conclusion 

We use data reconstruction techniques on the Voyager 2 data to study the radial evolution of 
the turbulent kinetic and normalized residual energy of the magnetic helicity and cross-helicity 
of the solar wind in “slow wind” conditions. We determine the probability density functions and 
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Figure 5: (a,b) Wave numbers spectra of the magnetic helicity density and of the magnetic energy /2 in Alfven 

units. These spectra are obtained by means of the Taylor hypothesis, see section]^ The mean speed is that computed 
over the whole time interval, see Table 1. For graphic clarity, in panel (a), for frequencies above 10“® Hz, not all the 
points of the spectra of the magnetic helicity are shown. The same in panel (b) above 10“^ Hz. (c,d) Relative normalized 
magnetic helicity CTmiK) = Kh,„lEm{K). The thick line indicates the cumulative relative helicity integrated from a given k 
to the highest observed wavenumber. 
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the spectra of the magnetic and cross-helicity at 5 and 29 AU from the Sun. We found a small 
level of Alfvenic fluctuations, which decrease with heliocentric distance, a possible indication 
of the evolution toward a turbulent state dominated by nonlinear interactions. The sign of the 
normalized cross-correlation remains positive, which indicates the slight prevalence of outward 
propagating perturbation waves. For the fluctuating field, the different time scales used in the 
statistical averages and cross-helicity spectra indicate that the outward and inward directions of 
propagation are equally probable in the whole range of scales we have analyzed. Together with 
the reduction of the normalized cross-helicity, this could suggest that, in the “slow wind”, solar 
wind turbulence is evolving toward a developed magnetohydrodynamic turbulence. However, 
data in 1989 show a sign change of polarization (the small scales are positively polarized) not 
present in 1979 data. Also, the polarization intensity is slightly higher in 1989. This result is not 
consistent with the hypothesized isotropization of the field predicted by models of homogeneous 
and isotropic magnetohydrodynamic turbulence ll^ . In future work we plan to extend this 
analysis to the outer part of the heliosphere and to the compressed solar beyond the termination 
shock. 
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